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For the first time the potential energy surfaces of actinide nuclei in the (20, 822, 830) deformation 
space are obtained from a multi-dimensional constrained covariant density functional theory. With 
this newly developed theory we are able to explore the importance of the triaxial and octupole shapes 
simultaneously along the whole fission path. It is found that besides the octupole deformation, the 
triaxiality also plays an important role upon the second fission barriers. The outer barrier as well as 
the inner barrier are lowered by the triaxial deformation compared with axially symmetric results. 
This lowering effect for the reflection asymmetric outer barrier is 0.5 ~ 1 MeV, accounting for 
10 ~ 20% of the barrier height. With the inclusion of the triaxial deformation, a good agreement 
with the data for the outer barriers of actinide nuclei is achieved. 


PACS numbers: 21.60.Jz, 24.75.+i, 25.85.-w, 27.90.+b 


Since the first interpretation of the nuclear fission by 
the barrier penetration [1], it has been a difficult task 
to describe this phenomenon theoretically. For example, 
in order to study the fission problem, one should first 
have very accurate information about the fission barrier; 
a 1 MeV difference in the fission barrier could result in 
several orders of magnitude difference in the fission half- 
life. Particularly, to explore the island of stability of 
superheavy nuclei (SHN) [2-4], it is more and more de- 
sirable to have accurate predictions of fission barriers of 
SHN [5-8]. 

Nowadays three types of models are used for calculat- 
ing fission barriers. During a long period the majority 
of these works is based on the macroscopic-microscopic 
(MM) models [5, 9-11]. The MM models make use of the 
Strutinsky shell correction method, allowing fast calcula- 
tions of multi-dimensional PES’s containing most of the 
important shape degrees of freedom. Until now it is still 
an important candidate for large scale fission barrier cal- 
culations based on the examination of multi-dimensional 
PES’s [10]. In recent years, the rapid development of the 
density functional theories (DFT) also makes it possible 
to calculate the fission barriers fully self-consistently [12- 
15]. There are mainly two reasons to apply DFT’s in the 
study of fission properties. First, many new functional 
forms and effective interactions are proposed with much 
better performances for the excited state as well as the 
ground state calculations [16-20]. Fission barrier calcula- 
tions are also helpful for developing these DFT’s. Second, 
in DFT’s much more shape degrees of freedom can be 
included self-consistently. For example, the symmetry- 
unrestricted Skyrme-Hartree-Fock-Bogoliubov model has 
been applied for the fission studies [21]. Besides these 
two types of models, there exist also methods intending 
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to combine the advantages of the MM and self-consistent 
models, such as the extended Thomas-Fermi method [22]. 


The double-humped fission barriers of actinide nuclei 
can be used to benchmark the predictive power of the- 
oretical models [13, 23-25]. Various shape degrees of 
freedom play important and different roles in the occur- 
rence and in determining the heights of the inner and 
outer barriers. For examples, it has long been known 
from MM model calculations that the inner fission bar- 
rier is usually lowered when the triaxial deformation is 
allowed, while for the outer barrier the reflection asym- 
metric (RA) shape is favored [26, 27]. Later on these 
points were also revealed in the non-relativistic [24] and 
relativistic [15, 28] density functional calculations, re- 
spectively. It is thus customary to consider only the tri- 
axial and reflection symmetric (RS) shapes for the in- 
ner barrier and axially symmetric and RA shapes for 
the outer one [5, 29, 30]. It has been pointed out that 
“there is no reason for a fissioning actinide nucleus not 
to penetrate all symmetry-breaking shapes on its way 
from the first (triaxial) to the second (mass-asymmetric) 
saddle” [31]. The non-axial octupole deformations are 
considered in both the MM models [32] and the non- 
relativistic Hartree-Fock theories [33]. However, a multi- 
dimensional structure of PES’s including both the triax- 
ial and RA shape degrees of freedom has not been ex- 
plored yet in the framework of covariant DFT. In this 
paper, we will investigate the influence of the triaxiality 
and the octupole shape on the PES’s all the way from 
the ground state to the fission configuration when both 
shape degrees of freedom are included simultaneously. To 
this end, not only as many as self-consistent symmetries 
should be broken, but also multi-dimensional constraints 
are needed [34]. 

To calculate the potential energy surfaces and fission 


barriers, in this work we use the covariant density func- 
tional theory (CDFT) [17, 19, 20, 35, 36]. By breaking 
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not only the axial [37, 38] but also the reflection sym- 
metries [39], we developed multi-dimensional constraint 
CDFT’s in which the functional can be one of the follow- 
ing four forms: the meson exchange or point-coupling 
nucleon interactions combined with the non-linear or 
density-dependent couplings [40]. If not specified, the 
functional form of the point-coupling nucleon interaction 
with non-linear self energy terms and the parameter set 
PC-PK1 [41] are used in this work. 

For the parametrization of the nuclear shape, we adopt 
the conventional ansatz in mean field calculations, 


At 


Pru = SAR 


(Qan) (1) 
where Q), is the mass multipole operators. When the 
axial and reflection symmetries are broken, the nuclear 
shape is invariant under the reversion of x and y axes. 
In other words, the intrinsic symmetry group is V4 and 
all shape degrees of freedom Sa, with u even, includ- 
ing the triaxial (u 4 0) and octupole (A = 3) deforma- 
tions, are possible. Irrespective with the self-consistent 
symmetries, the single particle wave functions and vari- 
ous densities are expanded on an axially deformed har- 
monic oscillator basis [38, 42]. In order to get a fast 
convergence of the results against the basis size, in the 
elongated direction more states are included in the ba- 
sis. Following Warda et al. [43], the basis is truncated as: 
nz/Q2+(2n,+|m])/Qi < Neut. Here nz, n1, and m are 
quantum numbers characterizing each state in the basis, 
Q- = MAX(1,b./bo), Q1 = MAX(1, b, /bo), and bo, bz, 
and b, are the oscillator lengthes. The calculated bind- 
ing energy of 74°Pu at B29 = 1.3 varies only about 130 
keV and 20 keV when Neut increases from 16 to 18 and 
from 18 to 20. This means a good convergence and such 
a truncation scheme with Nout = 16 ensures a 0.2 MeV 
accuracy for the deformation range we are interested in. 
In the present work, Neut = 16 (20) is used in the triax- 
ial (axial) calculations. More details of the convergence 
study will be given in Ref. [40]. The BCS approach is im- 
plemented in our model to take into account the pairing 
effect. Since it has been found that the BCS calculation 
with a constant pairing gap can not provide an adequate 
description of the fission barriers [14], we use a delta force 
for the pairing interaction with a smooth cutoff [41, 44]. 


We performed one- (1-d), two- (2-d), and three- 
dimensional (3-d) constrained calculations for the ac- 
tinide nucleus 74°Pu. In Fig. 1 we show the 1-d poten- 
tial energy curves (PEC) from an oblate shape with 820 
about —0.2 to the fission configuration with 829 beyond 
2.0 which are obtained from calculations with different 
self-consistent symmetries imposed: the axial (AS) or 
triaxial (TS) symmetries combined with reflection sym- 
metric (RS) or asymmetric cases. The importance of 
the triaxial deformation on the inner barrier and that of 
the octupole deformation on the outer barrier stressed 
by earlier studies [15, 24, 28] are clearly seen here: The 
triaxial deformation reduces the inner barrier height by 


E (MeV) 


FIG. 1. (Color online) Potential energy curves of °% Pu with 
various self-consistent symmetries imposed. The solid black 
curve represents the calculated fission path with V4 symmetry 
imposed, the red dashed curve that with axial symmetry (AS) 
imposed, the green dotted curve that with reflection symme- 
try (RS) imposed, the violet dot-dashed line that with both 
symmetries (AS & RS) imposed. The empirical inner (outer) 
barrier height Bemp is denoted by the grey square (circle). 
The energy is normalized with respect to the binding energy 
of the ground state. The parameter set used is PC-PK1. 


more than 2 MeV and results in a better agreement with 
the empirical datum; the RA shape is favored beyond 
the fission isomer and lowers very much the outer fis- 
sion barrier. Besides these features, we observe for the 
first time that the outer barrier is also considerably low- 
ered by about 1 MeV when the triaxial deformation is 
allowed. Again, a better reproduction of the empirical 
barrier height can be seen for the outer barrier. We note 
that this feature can only be found when the axial and 
reflection symmetries are simultaneously broken. 


How the PES of 74°Pu becomes unstable against the 
triaxial distortion can be seen much more clear in Fig. 2 
in which we show 2-d PES’s from calculations without 
and with the triaxial deformation. When the triaxial 
deformation is allowed, the binding energy of ?4°Pu as- 
sumes its lowest possible value at each (820, 830) point. 
At some points we get non-zero 22 values. That is, non- 
axial solutions are favored at these points than the axial 
ones. The triaxial deformation appears mainly in two 
regions in Fig. 2. One region starts from the first saddle 
point and extends roughly along the direction of the (39 
axis up to a very asymmetric shape with 639 ~ 1.0. In 
this region the values of 622 are about 0.06 ~ 0.12 cor- 
responding to y ~ 10°. The energy, especially the inner 
barrier height, is lowered by about 2 MeV. The other re- 
gion is around the outer barrier and the $22 values are 
about 0.02 ~ 0.03 corresponding to y ~2°. About 1 
MeV is gained for the binding energy at the second sad- 
dle point due to the triaxiality. In other regions, e.g., in 
the ground state and fission isomer valleys, only axially 
symmetric solutions are obtained. 

Next we examine the full 3-d PES of 74°Pu obtained 
from the newly developed multi-dimensional constraint 
CDFT. For simplicity, in Fig. 3 are shown only five typ- 
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FIG. 2. (Color online) Potential energy surfaces of 74°Pu in 
the (820, 830) plane from calculations (a) without and (b) with 
the triaxial deformation included. The energy is normalized 
with respect to the binding energy of the ground state. The 
numbers in (b) show the values of 622 at these points. The 
fission path is represented by a dash-dotted line. The ground 
state and fission isomer are denoted by full and open circles. 
The first and second saddle points are denoted by full and 
open triangles. The contour interval is 1 MeV. 


FIG. 3. (Color online) Sections of the three-dimensional PES 
of ° Pu in the (822, 830) plane calculated at 820= 0.3 (around 
the ground state), 0.6 (around the first saddle point), 0.9 
(around the fission isomer), 1.3 (around the second saddle 
point) and 1.6 (beyond the outer barrier), respectively. The 
energy is normalized with respect to the binding energy of the 
ground state. The contour interval is 0.5 MeV. Local minima 
are denoted by crosses. 


ical sections of the 3-d PES of 74°Pu in the (22, 830) 
plane calculated at 820= 0.3 (around the ground state), 
0.6 (around the first saddle point), 0.9 (around the fission 
isomer), 1.3 (around the second saddle point) and 1.6 (be- 
yond the outer barrier), respectively. Many conclusions 
can be drawn by examining these 3-d PES’s. First, the 
ground state and the fission isomer are both axially and 
reflection symmetric as what is shown in the 1-d PEC and 
the 2-d PES. But with the 3-d PES one can investigate 
the stability of 24°Pu against the 822 and (39 deforma- 
tions. One finds that the stiffness of the fission isomer is 
much larger than that of the ground state against both 
the 822 and $39 deformations. Second, while around the 
inner barrier the shape of 74°Pu is triaxial and reflec- 
tion symmetric, the second saddle point which is close to 
B29 = 1.3 appears as both triaxial and reflection asym- 
metric shape. Third, the triaxial distortion appears only 
on the top of the fission barriers. 

It has been pointed out that one may obtain spurious 
saddle points if only a small number of shape degrees 
of freedom are constrained, see, e.g., Ref. [5]. That is, 
the calculated fission path may jump from one valley 
to another and results in discontinuities in the lower- 
dimensional PES’s; in some cases, a continuous path 
may even cross a higher saddle point. Although the 
spurious saddle points may not be excluded completely, 
most of them can be avoided if (1) the obtained fission 
path keeps to be continuous in the energy as well as the 
most important shape degrees of freedom and (2) the 
results are examined by higher-dimensional calculations. 
We have carefully checked the full 3-d PES and found 
that the fission path enters and exits the triaxial con- 
figuration rather smoothly, which tells that no sudden 
jump is found and the 1-d (with the {29 deformation 
constrained and {22,339 deformations imposed) and 2- 
d (with (29,839 deformations constrained and the (22 
deformation imposed) calculations of the fission barriers 
may be well justified for 24°Pu. It is clear that the con- 
tinuity of the fission path found in a lower-dimensional 
constraint calculation is a necessary but not sufficient 
condition for locating the correct saddle point. In order 
to have a strictly definite conclusion, one certainly should 
carry out multi-dimensional constraint calculations with 
even higher-multipolarity deformations included. 

For the RS calculations, the triaxiality also lowers the 
fission path by a few MeV beyond the second saddle 
point. This point is illustrated by the dotted line in Fig. 1 
and the local minima with 339 = 0.0 in the G29 = 1.6 sub- 
figures of Fig. 3. However, it is relatively unimportant, 
because the RA fission is still the most favored one even 
when triaxiality is included. 

Guided by the features found in the 1-d, 2-d, and 3-d 
PES’s of ?4°Pu, the fission barrier heights are extracted 
for even-even actinide nuclei whose empirical values are 
recommended in RIPL-3 (see Table XI in [45]). The em- 
phasis is put on the influence of the triaxial deformation 
on the two fission barriers. 

As it has been shown previously, around the inner bar- 
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FIG. 4. (Color online) The inner (B}) and outer (Bẹ) barrier 
heights of even-even actinide nuclei. The axial (triaxial) re- 
sults are denoted by open (full) symbols. The empirical values 
are taken from Ref. [45] and represented by grey squares. 


rier an actinide nucleus assumes triaxial and reflection 
symmetric shapes. Thus in order to obtain the inner fis- 
sion barrier height we can safely make a one-dimensional 
constraint calculation with the triaxial deformation al- 
lowed and the reflection symmetry imposed. In Fig. 4(a) 
we present the calculated inner barrier heights Bi and 
compare them with the empirical values. It is seen that 
the triaxiality lowers the inner barrier heights of these 
actinide nuclei by 1 ~ 4 MeV as what has been shown 
in Ref. [15]. In general the agreement of our calculation 
results with the empirical ones is very good with excep- 
tions in the two thorium isotopes and ?°°U. For 72°Th and 
232Th, the calculated inner barrier heights are smaller by 
about 2 or 1 MeV than the empirical values depending 
on whether the triaxial deformation is allowed or not. In 
these two nuclei, the outer barrier is higher than the inner 
one. This may result in some uncertainties when deter- 
mining empirically the height of the inner barrier which 
is not the primary one [12]. Similar results for 7°°Th 
and ?3?Th were obtained from the Skyrme-Hartree-Fock- 
Bogoliubov model [12] and very small inner barrier height 
was got for ?3?Th in Ref. [15]. For ?°8U, Bi from the ax- 
ial calculation agrees with the empirical value very well. 
The triaxiality reduces the barrier height by about 1.5 
MeV, thus bringing a discrepancy which was similar to 
the result in Ref. [15]. 


To obtain the outer fission barrier height BẸ, the sit- 
uation becomes more complicated because more shape 
degrees of freedom have important influences around the 
outer fission barrier. For example, the inclusion of the 
reflection asymmetric shape makes it possible to have in 
the (20, 830) plane two or more competing fission paths 
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FIG. 5. (Color online) (a) One-dimensional potential energy 
curve E ~ 629 and (b) two-dimensional potential energy sur- 
face E ~ (620, B30) of °“8Cm the outer barrier region with the 
axial symmetry imposed in the calculation. In both figures, 
the energy is normalized with respect to the binding energy 
of the ground state. The fission path I/II is represented by 
full/dash dot lines and the corresponding saddle point is de- 
noted by up/down triangles. The fission isomer is denoted 
by open circle. In (a), the empirical outer barrier height is 
depicted by dash-dotted line. In (b), the contour interval is 
0.5 MeV. 


with different octupole deformations. In consequence, 
one often observes in the 1-d PEC two or more fission 
paths. This happens in 244?46Pu and ?44246.248Cm in 
the present study and we present a typical example in 
Fig. 5 for 24°Cm. In this figure, one finds that there are 
two fission paths both in the 1-d Æ ~ 69 curve and in 
the 2-d E ~ (629, 639) PES. One path denoted by “I” 
favors shapes with larger octupole deformations and the 
other denoted by “II” favors less RA shapes. In such 
nuclei, it is not safe to perform a 1-d constraint calcula- 
tion in order to get BP. Thus we first assume the axial 
symmetry and make a 2-d calculation in the (20, 830) 
plane from which we can approximately identify the low- 
est fission path 619"°st (G29) and the location of the sec- 
ond saddle point. Then along this fission path, we per- 
form a 1-d (9-constraint calculation with the triaxial 
and octupole deformations allowed. At each point with 
B20, the initial deformations are taken as pił- = 0 and 

ini. = Blowest (Bo). In this 1-d PEC, we can locate the 
second saddle point and extract the outer barrier height 
for each nucleus. 

In the lower panel of Fig. 4 we show the results of outer 
barrier heights BP and compare them with empirical val- 
ues. For most of the nuclei investigated here, the triaxi- 
ality lowers the outer barrier by 0.5 ~ 1 MeV, accounting 
for about 10 ~ 20% of the barrier height. One finds that 


201708.00254v1 


chinaXiv 


our calculation with the triaxiality agrees well with the 
empirical values and the only exception is 24°Cm. From 
the calculation with the axial symmetry imposed, the 
outer barrier height of ?48Cm is already smaller than the 
empirical value. The reason for this discrepancy may be 
related to that there are two possible fission paths beyond 
the first barrier, as seen in Fig. 5. For the path I with 
a lower saddle point from which we get the outer fission 
barrier height, the barrier is very wide and for the path I 
with a higher saddle point, the barrier is relatively nar- 
row. Therefore the empirical value of the outer fission 
barrier height may not be easily extracted for the follow- 
ing two reasons: (i) There must be a strong competition 
between the two fission paths; (ii) When the empirical 
value of the outer barrier height is evaluated, it is usually 
assumed that the second barrier is in an anti-parabolic 
shape with a fixed and smaller width [45]. 


We also examined the parameter dependency of our 
results. The lowering effect of the triaxiality on the outer 
fission barrier is also observed when parameter sets other 
than PC-PK1 are used. 


In summary, a multi-dimensional constrained covariant 
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